Theory of quantum noise detectors based on resonant tunneling 
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We propose to use the phenomenon of resonant tunneling for the detection of noise. The main idea 
of this method relies on the effect of homogeneous broadening of the resonant tunneling peak induced 
by the emission and absorption of collective charge excitations in the measurement circuit. In 
thermal equilibrium, the signal to noise ratio of the detector as a function of the detector bandwidth 
(the detector function) is given by the universal hyperbolic tangent, which is the statement of the 
fluctuation-dissipation theorem. The universality breaks down if non-equilibrium processes take 
place in the measurement circuit. We propose the theory of this phenomenon and make predictions 
for the detector function in case when non-equilibrium noise is created by a mesoscopic conductor. 
We investigate measurement circuit effects and prove the universality of the classical noise detection. 
Finally, we evaluate the contribution of the third cumulant of current and make suggestions of how 
it can be measured. 

PACS numbers: 72.70.+m, 42.50.Lc, 74.50.+r, 73.23.Hk 
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I. INTRODUCTION 

Universalities play an important role in physics, be- 
cause they point to fundamental laws and properties, 
such as symmetries, topology, scaling behavior, and oth- 
ers. Moreover, when broken, they open a door to new 
physics. Here we wish to consider one example that is im- 
portant in the context of the present work. Recently, fol- 
lowing the suggestion of Kane and Fisher, 1 experiments 
on shot noise in quantum Hall systems 2 ' 3 directly mea- 
sured fractional charge of Laughlin quasiparticles. The 
interpretation of these experiments invokes a simple ar- 
gument that weak quasiparticle tunneling is an uncorre- 
cted Poisson process which is described by the Schottky 
formula S — q(I), where (J) is the average tunneling cur- 
rent, S is the zero- frequency noise power of the tunneling 
current, and q is the fractional charge of quasiparticles. 

More rigorously, the Schottky formula follows from the 
fluctuation dissipation theorem (FDT), which states that 
when a tunnel junction weakly connects two metallic 
reservoirs, the following relation generally holds 4,5 



v+sv 



q(I)/S = tanh(A/V2fc s T), 



(1) 



where A/x is the electro-chemical potential difference ap- 
plied to the barrier. This relation is a generalization of 
the well-known Callen-Welton FDT which connects the 
noise power and the linear response coefficient, 8 and fol- 
lows from the argument similar to the one used in the 
linear response theory. This implies the universality of 
the relation (1), i.e. it holds independently of the char- 
acter of the interaction, spectrum of quasiparticles, the 
geometry of a tunnel junction, etc. It is easy to see that 
the Callen-Welton theorem and the Schottky formula are 
the two limits of the relation (1). 

Here we present a simplified derivation of (1), based 
on the "Golden rule". 9 Quantum mechanical transition 
rate between to energy states E n and E m is given by 
W mn = 2%5(E n - E m )\A mn \ 2 , where A mn = (E m \A\E n ) 
is the matrix element of the tunneling amplitude A. 
Then the average current can be evaluated as (I) = 
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FIG. 1: The measurement circuit contains a mesoscopic sys- 
tem which creates non-equilibrium noise and has the conduc- 
tance Gm, the detector consisting of a tunnel junction shunted 
by the load resistor, Gl, and the capacitor C. The voltage 
bias Vm and Vd is applied to the system and to the tunnel 
junction, respectively. Fluctuating current through a meso- 
scopic system is accumulated on the capacitor and creates the 
fluctuating potential SV across the tunneling barrier. 



lUmn W mn (p n -p m ), where p n = (E n \p\E n ) is the diag- 
onal matrix element of the density operator, i.e. the prob- 
ability to find the system in the state E n . When tunnel- 
ing is weak, forward and backward tunneling transitions 
are independent Poisson processes with the dispersions 
of fluctuations equal to mean currents. Therefore, the 
total noise power is equal to S — q 2 J2 mn W mn (fi n + Pm)- 
In equilibrium p n cx exp(— S n /fcflT), so that p n = p m , 
and the current vanishes. If the potential difference A/i 
is applied between the leads which are locally at equi- 
librium, then the density matrix acquires the grand- 
canonical form p(A/j) = p(0) exp(A/iiV//csr), where 
N is number of electrons in one of the leads. Since 
the tunneling amplitude changes the number of parti- 
cles in this lead by one, then obviously one can write 
p m = exp(A/i//csr)p„, which immediately gives the re- 
lation (1). 

From the derivation of the FDT it is obvious that non- 
equilibrium processes in reservoirs play a special role, 
since they may lead to a deviation from the universal 
relation (1). The goal of the present paper is to investi- 
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gate this phenomenon in the case when noncquilibrium 
processes take place in the electrical (measurement) cir- 
cuit to which the tunnel junction is attached. In Fig. 1 we 
draw its simplified version that contains essential parts: 
The mesoscopic system which creates non-equilibrium 
noise and has the conductance Gm, the detector con- 
sisting of a tunnel junction shunted by the load resistor, 
Gl, and the capacitor C. One of the important results 
of our paper is that the FDT breaks down in a minimal 
way, so that some properties of the current-to-noise ra- 
tio, which contain an important information about non- 
equilibrium processes in the leads, retain their univer- 
sality. This leads to the idea to use tunnel junctions as 
on-chip detectors of non-equilibrium noise, which we in- 
vestigate below in details. 

The measurement circuit has been proven to play an 
important role in the physics of the noise detection with 
the standard measurement technique 10-12 and with the 
help of on-chip noise detectors. 13 It has been established 
that in the long time (Markovian) limit the backaction 
of the measurement circuit on the system leads to "cas- 
cade corrections" to statistics of noise. 14,15 In order to 
quantify the circuit effects one solves the Kirhgoff (the 
current conservation) law for the fluctuations of the cur- 
rent through the mesoscopic system 51 m and through the 
load resistor 51 d , and the voltage fluctuations 5V on the 
capacitor: 

SV(lj) = Z(u)[5I M (u) + 5I L (u)]. (2) 

The circuit impedance is given by 

Z(w) = R/(l - iujtc), (3) 

where R = 1/(Gm + Gl) is the differential circuit resis- 
tance, and tc = RC is the circuit response time. The 
Eq. (2) describes the effect of the system current fluctua- 
tions via the circuit on the tunnel junction which directly 
detects potential fluctuations. The normalized circuit re- 
sistance TZ — G a R, where G = e 2 /2ir is the conductance 
quantum, plays a role of the dimcnsionlcss coupling con- 
stant, which parametrizes the strength of the circuit ef- 
fects. In the present work we assume that coupling is 
weak, R< 1. 

Quantum noise detectors, the main operating principle 
of which is based on the resonant tunneling in a two-level 
system, where investigated experimentally and theoreti- 
cally in a number of previous works. 16-21 Here we con- 
sider two different detectors of this type. The first one, 
the double-dot (DD) detector of quantum noise theoreti- 
cally analyzed by Aguado and Kouwenhoven, 17 is shown 
in Fig. 2. It consists of two quantum dots, which are 
strongly coupled to leads and weakly coupled to each 
other. To the lowest order in inter-dot coupling the elec- 
tron transport in the DD detector occurs via inelastic 
transitions between nearest energy levels of two dots. 
These transitions are assisted by the emission (absorp- 
tion) of the energy e to (from) the circuit, where e is the 
inter-dot level distance. Away from the resonance the 



average current through the DD detector is given by 17 
(/) ~ Sm{s)/£ 2 , where Sm(e) is the spectral density of 
the non-symmetrized correlator of the system current. It 
is easy to see that the parameter e plays a role of the 
bandwidth of the detector. 

The operating principle of the second detector, based 
on the telegraph process (TP detector, see Fig. 3), is 
slightly different. It contains two weakly coupled quan- 
tum dots, which are electrically isolated from the circuit, 
but capacitively coupled to it. Fluctuations of the po- 
tential on the capacitor, caused by the current noise in 
the mesoscopic system, lead to rare electron transitions 
between two dots which change the electrical charge of, 
say, the left dot randomly in time. When left dot in- 
teracts with a nearby quantum point contact (QPC), it 
randomly switches the QPC current between two levels, 
Id and I u , leading to the telegraph process. The average 
QPC (detector) current (Id) is a monotonic function of 
the inter-dot level distance e, which changes from one 
current level to the other, depending on the average oc- 
cupation of the left dot. Thus the QPC acts as a sensi- 
tive electrometer of the occupation of the quantum dot, 
the principle demonstrated in early work [22] and subse- 
quently elaborated in recent experiments, where the real 
time detection of single-electron tunneling, 23-28 the mea- 
surement of counting statistics, 29-31 and the information 
backaction of a detector 32 have also been shown. 

We denote with D(e, A) the current-to-noise ratio for 
the DD detector and call it the detector function. Al- 
though D depends on the bias voltage Vd, we choose to 
represent it as a function of the energy A = eVu — £ 
of the electron-hole pair created in the leads by the ele- 
mentary tunneling event. This energy parametrizes the 
asymmetry of the detector, because in case A = 0, or 
equivalently e = cVd, there is no difference between left 
and right dot of the DD system. Below we prove an im- 
portant fact that the average current through the QPC 
of the TP detector, after a proper normalization [see Eqs. 
(26) and (27)], is given by the symmetric detector func- 
tion D(e, 0). When the circuit is at thermal equilibrium, 
D(e,0) = tanh(e/2fc B r) according to the FDT. 

The physics of quantum noise detection is quite rich 
thanks to a number of energy scales that determine the 
dynamics of entire system. While these energy scales are 
not important in the case of equilibrium circuit, because 
the FDT holds and leads to the universality, they start 
to play an important role away from equilibrium. First 
of all, it is an effective temperature of the noise source, 
0, which is formally defined by Eq. (18). It has a mean- 
ing of the energy provided by the system and the load. 
Alternatively, one can think of the correlation time 1/17 
of the noise source. Second important parameter is the 
detector bandwidth e, introduced earlier. Third, the cir- 
cuit itself is characterized by the response time tc and 
corresponding energy scale 1/tc (we set h = 1). Finally, 
the asymmetry of the DD detector is characterized by 
the energy A. 

In the weak coupling limit, which we consider through- 



3 



out the paper, where the dimensionless circuit impedance 
1Z is small, the detector operates at the Gaussian point, 
i.e. the contribution of high order cumulants (irreducible 
moments of noise) is small. The physical reason for this 
is that in the limit 1Z <C 1 the detector only weakly inter- 
acts with the the noise source, therefore it has to operate 
for a relatively long time interval of the order of 1/1ZQ. 
in order to accumulate a sufficient information about the 
noise. During this time interval many fluctuations con- 
tribute to the detector signal, so that by virtue of the cen- 
tral limit theorem the resulting noise becomes Gaussian. 
In Sec. VII we show how the third cumulant, which is the 
simplest characteristics of the non-Gaussianity, neverthe- 
less can be extracted from the detector output signal. 

The new energy scale Tq arises in the weak coupling 
limit due to the effect of homogeneous level broadening: 
Close to the resonance, e — ► 0, the interaction of the 
detector with the circuit becomes effectively strong, and 
inelastic transitions in the detector are assisted by mul- 
tiple photon absorption and emission processes. As a 
result, the detector signal at this point acquires a peak 
as a function of e of the width Tq -C ^- The shape of the 
peak depends on the circuit details. We distinguish two 
limiting cases, depending on the circuit response time 
tc- In the "fast" circuit limit, IZQtc -C 1, the peak 
has a Lorentzian shape and the width Tq — 2n1Z£l, see 
Eq. (21). In the "slow" circuit limit, TZflrc > 1, the 
peak acquires the Gaussian shape (23) with the width 
Tq = where Ec = e 2 /C is the Coulomb charg- 

ing energy of the circuit. 

Depending on energy scales following regimes can be 
distinguished. In the quantum noise detection regime, 
e <~ Q, the detector signal is due to the inelastic tunnel- 
ing with the absorption or emission of a single photon of 
the energy e. The probability of this process is given by 
Eq. (16). In case when the circuit is driven away from 
equilibrium by a coherent mesoscopic conductor the sym- 
metric detector function is given by Eq. (31). For a low 
circuit impedance, GmR <C 1, this expression simplifies 
and we obtain the result (32). The results are summa- 
rized in Fig. 4. In the case of strongly asymmetric detec- 
tor, A 3> f2, the detector function takes the equilibrium 
tangent form (33), which is however shifted by the energy 
Em given by Eq. (34), which can be viewed as the noise 
rectification effect. 

In the classical noise detection regime, e <C fi, the de- 
tector function is simply linear in e [see Eq. (35)], with 
the slope determined by the effective noise temperature 
Q.. Thanks to this universality, there is no need to spec- 
ify the mesoscopic system which is measured. Close to 
the resonance, e ~ Tq, the inelastic tunneling becomes 
non-perturbative despite the small parameter 1Z, and the 
P(E) function acquires a peak of the width Tq. The 
shape of the peak depends on the circuit details, see Eqs. 
(21) and (23). Nevertheless, the detector function retains 
its universal form (35), so it can be used to extract the 
noise temperature. 

We evaluate the small contribution of the third cumu- 



lant of the system current in the classical (Markovian) 
limit and find that it slightly shifts the zero of the detec- 
tor function (43) by the energy E3 which is proportional 
to the third cumulant. The coefficient depends on the 
circuit response time tc and is evaluated in the case of 
fast and slow circuit, see Eqs. (46) and (49). The total 
third cumulant of the system current contains cascade 
corrections, which depend on the circuit response time. 
In the case of fast circuit the cascade corrections are given 
by Eq. (45), i.e. they are those introduced by Nagaev in 
Ref. [14]. In the slow circuit case the detector measures 
equal time fluctuations of the potential on the capacitor, 
and the cascade corrections in this case are given by Eq. 
(48), as predicted in Ref. [13] and measured in Ref. [12]. 
We finally note, that the third cumulant of current may 
be extracted from the shift of the detector function us- 
ing the technique recently introduced in experiments on 
the mesoscopic threshold detectors. 33,34 The universality 
of the detector function in the classical noise detection 
regime, proven in Sec. VI, may become crucial for the 
success of this procedure. 

The rest of the paper is organized as follows. After re- 
viewing the P(.E)-theory of tunneling in Sec. II, we focus 
on the Gaussian noise case in Sec. Ill and classify the 
measurement circuit effects according to the circuit re- 
sponse time. In Sec. IV we analyze quantum noise detec- 
tors based on the resonant tunneling effect and connect 
the detector function D(e, A) to the current-to-noise re- 
lation for tunnel junctions. We use results of the P(E)- 
theory in Sec. V to calculate the detector function in 
quantum and classical noise detection regimes. In Sec. 
VI we prove that the detector function is universal in 
the classical noise detection regime, i.e. it is independent 
of the measurement circuit details. Finally, in Sec. VII 
we investigate the third cumulant contribution to the de- 
tector function including the circuit cascade corrections. 
The section VIII outlines further directions of research. 

II. REMINDER ON P(£)-THEORY OF 
TUNNELING 

The purpose of this section is to remind essential steps 
of the P(P)-theory of photon-assisted tunneling. 35 In ad- 
ditional, we extend the theory in order to take into ac- 
count weak non-Gaussian effects in noise. The tunnel 
junction, attached to two metallic leads, is described by 
the Hamiltonian: 

H = J2 E k(clc k + dld k ) + H T , (4) 
k 

where Cfe and dk are the electron operators in the left and 
right lead, respectively, and Ht is the tunneling Hamil- 
tonian. It can be written as 35 

H T = A + A\ A = e^Y. T p* d l c ^ ( 5 ) 

pk 

where the amplitude A transfers the electron from left to 
right, and the phase factor e 1 ^ changes the charge on the 
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capacitor by — e. The last fact follows from the charge 
quantization [0, Q] = ei. Then the charge Hamiltonian 
He = Q 2 /2C generates the equation of motion for the 
phase operator: <fi = eSV. We thus assume that the inter- 
action of electrons with the collective charge excitations 
in the electrical circuit is generated solely via tunneling. 

Next we evaluate the average tunneling current 
(I D ) and the zero- frequency noise power Sp = 
J dt(5In(t)5lE)(0)) . We define the tunneling current op- 
erator as Id = edN^/dt = ie(A — A'), where Nl — 
J2k c l c fe i s tne number of electrons in the left lead. To 



leading order in the tunneling Hamiltonian (5), we can 
write 

(Id) = e J dt([A(t),AH0)]}, (6a) 
S D = e 2 J dt({A(t),A^(0)}). (6b) 

Substituting A from Eq. (5) to Eqs. (6), and tracing out 
electronic operators, we finally obtain: 



(Id) = 2neJJ dEdE 1 v R {E)v L {E'){P LR (E - E' + eV D )f(l - /') - P RL (E' -E- eV D )f'(l - /)}, (7a) 
dEdE'v R (E)v L (E'){P RL (E -E' + eV D )f(l - /') + P LR (E' -E- eV D )f'(l - /)}, (7b) 



S D = 2ne' 



where / = /f(E) and /' = /f(-E') are the equilibrium 
distributions in the leads, and vl and v R are the elec- 
tronic densities of states. Here Plr(E) and P R l(E) are 
the probability distributions of the emission (absorption) 
of a collective charge excitation of energy E, caused by 
inelastic tunneling of an electron from right to left lead 
and vice versa: 



Plr(E) 



1 

2^ 



Prl(E) = ±J dte lEt (e~ l 



dte iEt( e it(t) e -i<KO)) t ( 8a ) 

(8b) 



0(t)g#(O) 



In general, the phase correlation functions in (8) can be 
expanded in terms of the noise cumulants. However, ev- 
ery cumulant comes with extra coupling constant 1Z <SC 1. 
Therefore, we keep only first two nonvanishing cumulants 



J 2 (t) = i(^ 2 (t)-20(t)0(O) + 2 (O)) 

Mt) = \(^{t)-?,<p 2 {t)m 



and write 



Plr(E) 
Prl(E) 



+30(t)0 2 (O)-^(O)) 

— f dte iEt ~ j2( >' t) ~ ij3< > t \ 
2tt J 

— ( dte lEt ^ j2{t)+t - h{t) 
2ttJ 



(9) 
(10) 

(11a) 
(lib) 



We postpone the discussion of the third cumulant effect 
till Sec. VII and for a moment assume that the noise is 
Gaussian. 



III. GAUSSIAN NOISE 

We now set J3 = and write Plr = Prl = P, where 

P(E) = — [ dte iEt - J 'W, (12) 
2tt J 



Note that in Eq. (9) each term of the form (<fi 2 ) con- 
tain a classical contribution, which in the long time limit 
is proportional to time. 36 This is a consequence of the 
Brownian motion of the phase "pushed" by a fluctuating 
potential. However, these potentially dangerous terms 
cancel, and Eq. (9) can be rewritten in the form 

J 2 = i((A^) 2 ) + i([0,A0]), A0 = m - 0(0), (13) 

so that it does not contain divergences. In this equation 
the first term can be interpreted as a classical contribu- 
tion, which is proportional to time in the long time limit, 
and the second term is pure quantum. Using Eq. (2), we 
obtain 



Mt) = Go 



duo S(lo) J 
uj 2 + rj 2 



Z( W )| 2 (l-e— ^0 (14) 



where S(uj) = Sm(w) + Sl{u) is the power of the total 
noise created in the circuit, and Sm and Sl are the non- 
symmetrized correlators of the mesoscopic system and of 
the load resistor 



S M (u) = J dte iut (6I M (t)SI M (0)), 
SlH = I dte wt (5I L {t)5I L {0)), 



(15a) 
(15b) 



Next we note that in the weak interaction case 1Z = 
Go-R <SC 1 considered here, J 2 (t) is usually small. For 
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instance, in equilibrium RS(ui) = 2fcsT, so that for 
t ~ 1/ksT the correlator given by Eq. (14) can be 
roughly estimated as J 2 ~ 1Z. Therefore, we expand 
the exponential on the right hand side of Eq. (12) and 
obtain: 



P(E) = P Q S(E) + G \Z(E)\ 2 S(E)/E 2 



(16) 



where Po is the probability of the elastic process, fixed by 
the normalization J dEP(E) = 1. The probability of the 
inelastic process is proportional to the non-symmetrized 
correlator S (E) , 17 ' 37 ~ 39 and at relatively large energies is 
sensitive to quantum fluctuations. 

Special care however has to be taken about the long 
time limit in Eq. (14), since growing with time classical 
contribution to J 2 may compensate smallness of 1Z. The 
Fourier integral cuts off a small region around u> = 0, 
where the noise is classical, and the noise power can be 
approximately replaced with 5(0). The important note is 
in order: Quantum effects, which lead to the interaction 
induced suppression of tunneling, i.e. to so called dynami- 
cal Coulomb blockade effect, 35,40 are not neglected. They 
are fully taken into account in Eq. (16), and subsequently, 
in Sec. V. However, at the energy scale of interest here 
their contribution to the long time asymptotic is small. 
We now focus on the long time limit and consider the 
cases of fast and slow circuit, depending on the circuit 
response time tc- 



A. Fast circuit 

We first assume that the relevant time scale is longer 
than tc, and therefore set Z(u) = R. From Eq. (14) we 
find: 

J 2 (t) = 2n1in{\t\ + id ul S(0)/S(0)sign(t)}, (17) 
where the energy scale fi is the circuit noise temperature: 



n = (1/2)RS(0), 



(18) 



Note that although the interaction is weak, 1Z <C 1, in 
the long time limit \t\ ~ 1/1ZCI the exponential in Eq. 
(8) cannot be expanded. We then use the result (17) and 
obtain: 



P(E) 



2nn 



E 2 + (2kKQ)< 



[1 + £790,5(0)75(0)], (19) 



which is consistent with the result (16) in the limit E ^> 

nn. 

Thus we find that in the limit, \t\fl ~ 1/11 > 1, the 
multiple photon processes lead to the broadening of the 
(5-function in Eq. (16), so that it is replaced with the 
Lorentzian peak. One can now use Eq. (3) to check that 
the assumption Z(uj) = R is justified if IZQtc <C 1. This 
means that the response of the circuit to current fluctu- 
ations is instantaneous, and the phase fluctuations are 
Markovian on the time scale of interest. 



The asymmetry of P(E) given by the second term 
in Eq. (19) is weak: Ed u S/S <~ TZ. Interestingly, 
the expression d u S(0) = J dt(it/2)([I(t), 1(0)]) coincides 
with the Kubo formula for the differential conductance, 
l/R = d v (I). Therefore we obtain d u S(0) = 1/R, and 
alternatively, 



du,S(0) 
5(0) 



1 

2ft 



(20) 



Thus one can express the asymmetry in (19) in terms of 
the noise temperature alone: 



P(E) 



K(2n + E) 

e 2 + (2Trnny 



nriTc < 1. 



(21) 



B. Slow circuit 

Next we consider the opposite limit, 1Z£It c » 1, when 
circuit responds slowly to current fluctuations. In this 
case it is the singularity in Z(u) cuts off the integral in 
(14) at small frequencies uo ~ 1/rc- Using the impedance 
(3) and the relations (18) and (20), we obtain: 



J 2 (t) = ir1Z(n/T C )t 2 + iir1l(l/Tc)t 



(22) 



The first term in this equation has a simple interpre- 
tation. We note that it can be also obtained by con- 
sidering the phase (f> as classical variable, and writing 
cj>(t) — (f)(0) = eSVt, because the variation of the poten- 
tial is slow. Then, the equation (13) leads to J 2 (t) = 
(e 2 /2)((SV) 2 )t 2 , which together with Eq. (3) gives the 
first term in (22). Thus the phase correlator is deter- 
mined by the equal time correlator of the potential. We 
will rely on this interpretation in Sec. VII. 

The second term in the equation (22) has a quantum 
nature. It slightly shifts the energy in P(E) given by 
Eq. (12) and leads to the asymmetry of the distribution. 
This term is small, therefore the Fourier transform in (8) 
can be written as 

where Ec = e 2 /2C. Thus we see that the dissipative 
properties of the circuit, determined by the resistance R, 
do not enter the final result. This is related to the fact, 
that slow fluctuations of charge on the capacitor obey the 
low of the equipartition of energy, (1/2C)((<5Q) 2 ) = fi/2, 
as it were in equilibrium. 

We remark that in the intermediate regime IZflrc ~ 1 
the exact shape of the zero energy peak in P(E) is more 
complex and depends on details of the circuit. Neverthe- 
less, as we show in Sec. VI, the fluctuation-dissipation 
relations remain insensitive to these details. Finally, we 
also note that these our results on the asymmetry in 
P(E) were published in Ref. [41]. Recently the asym- 
metry was found in the experiment [42] and theoretically 
discussed in Ref. [43]. 
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FIG. 2: The DD detector operates as shown on the upper 
panel. The absorption of the quantum of the collective charge 
excitation of the circuit leads to the inelastic electron transi- 
tion between two weakly connected quantum dots. Because 
dots are strongly connected to two metallic reservoirs, multi- 
ple random transitions generate current through the detector 
and the current noise. Lower panel show the energy diagram 
of the detector and most important parameters. 

IV. QUANTUM NOISE DETECTORS 




FIG. 3: The TP detector consist of a double-dot system, 
which is capacitively coupled to a QPC. When one of the 
dots is charged, it pinches off the QPC and thus changes the 
current through it from the upper level I u to the lower level 
Id- Random interdot transitions, caused by the emission and 
absorption of the collective charge excitations of the circuit, 
lead to random switching of the QPC current. The result- 
ing telegraph process is shown on the upper right panel. The 
average current through the QPC, shown on the lower right 
panel, develops a smooth crossover between two current levels 
as a function of the DD level distance e. 



We have briefly discussed two types of quantum noise 
detectors in the introduction. Here we analyze them in 
details and show that their properties are determined by 
the P{E) function, obtained in previous section. Starting 
with the DD detector, we first assume that tunneling 
between two dots is a weakest process. In this simple 
case the transport can be described by lowest order in 
tunneling, so that the result (7) of previous section fully 
applies. Moreover, a weak coupling of the dots to the 
reservoirs leads to the broadening of the dot levels, so 
that the densities of states acquire a Breit-Wigner form 
v a = {T a /n)/[(E - E a f + Til a = L,R. 

If the noise temperature is small, so that To < T a , 
then the elastic transport dominates the photon-assisted 
inelastic transitions. In this case the left and right lead 
are approximately at thermal equilibrium, and the FDT 
holds. The most efficient noise detection takes place for 
a relatively strong noise in the circuit, Tji > r a , when 
the homogeneous level broadening dominates the quan- 
tum effect. In this case Breit-Wigner resonances can be 
replaced by delta functions, v a = 5(E — E a ), a = L,R, 
where El and Er are the energies of dot levels counted 
from the local Fermi level in the left and right lead, see 
Fig. 2. Substituting delta functions to Eqs. (7) and us- 
ing /'(l - /) = /(l - /') e ( £ «-^)/ fc B T for the current to 
noise ratio we obtain the following function: 



e -M^D(s,A)- P{£)eA/kBT 
Sd 



P(-e) 



P( e ) e A/fc B T + P (_ £ )- 



(24) 



Here the tunable level distance e = Er — El + eVo is 
the detector bandwidth, and the energy of the electron- 
hole pair A = El — En parametrizes the asymmetry 
of the detector. The detector function D(e, A) will be 
analyzed in details in the next section. Below we show 
that the properties of the TP detector are determined by 
the symmetric variant of this function, -D(e, 0). 

We evaluate the average current through the QPC ca- 
pacitively coupled to the DD, see Fig. 3. Switching of 



the DD from one state to another changes the current 
Id through the QPC from the low level Id to the high 
level /„. The probability to find the DD in the lower and 
upper state are given by P d = j d / (j u + j d ) and P u = 
lu/(lu + Id), where id and i M are the switching rates. 
Then the average current is given by (Id) = I U P U + IdPd 
and is equal to 



{Id) 



lulu + hid 

In + Id 



(25) 



It is convenient to rewrite the detector current (Id) m 
the dimensionlcss form: 



_ 2(I D ) - (Iu+Id) _ lu-ld 



Id 



lu +ld 



(26) 



so that it acquires the maximum value Td — 1 when 
lu S> id and the upper level is occupied, P u = 1, and 
Td = — 1 in the opposite case, when mostly the lower 
level is occupied. 

Next we assume that one of the dots is strongly cou- 
pled to the circuit capacitor. Then with a good approx- 
imation switching of the DD changes the charge of the 
capacitor by the value e, so that the interdot coupling is 
proportional to e l< ^. 44 Assuming the interdot coupling is 
weak compared to the width of levels, one can evaluate 
the switching rate using the Golden rule approximation 
with the result j u oc P(s) and id oc P(—e), where e is 
the DD level distance. Therefore, using the result (26) 
we obtain: 



1d(s) = 



P(e) - P(-e) 
P(e) + P(-e) 



D(e,0), 



(27) 



i.e. the normalized average current through the QPC as 
a function of the tunable level distance is given by the 
symmetric variant of the detector function. 
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V. FLUCTUATION-DISSIPATION RELATIONS 

In this section we investigate how non-equilibrium pro- 
cesses in the circuit lead to a breakdown of the FDT. We 
first focus on the inelastic regime e > IZtt, where Eq. 
(16) applies, and later consider the classical regime de- 
scribed by Eqs. (21) and (23). Substituting Eq. (16) to 
the definition (24), we obtain 



S(e)e 



A/k B T 



S(- 



(28) 



where, we remind, e = eVo — A is the interdot level 
distance. Thus all the circuit details cancel from the 
final result, and the exact form of the function D is de- 
termined solely by a non-symmetrized correlator of the 
current fluctuations in the circuit. 

In order to make further progress, we have to specify 
the model of the current source. The load resistor may 
be considered as a macroscopic system that creates an 
equilibrium current noise. Non-equilibrium processes are 
generated by the mesoscopic system alone. An interest- 
ing and experimentally important example of the meso- 
scopic system is a coherent mesoscopic conductor, which 
is fully characterized by a set of transmission eigenval- 
ues, T n , n = 1, . . . , N. Using the scattering theory, 45 one 
obtains the following expression for the non-symmetrized 
current correlator 

S m (oj) = GoY, {^F(uj)+T n (l-T n ) 

x [F(u + eV M )+F(w-eV M )]\, (29) 



where 



F(w) 



1 



-u/k B T ' 



(30) 



and we assumed that transmission eigenvalues T n are en- 
ergy independent. Using F(u>) — F(—u>) — 2to we now 
check that indeed Sm(^>) — Sm{—^) = ^ljGmi where the 
conductance Gm = Go~Y^ n T n . The same relation obvi- 
ously holds for the macroscopic resistor. 

In equilibrium Eq. (29) gives S = S M + S L = 2(G M + 
Gl)F(uj), a well known result for the non-symmetrized 
noise power. It satisfies the detailed balance relation 



S(-co) 



j/ki 



t S(lu). Substituting this relation to 



Eq. (28), we arrive at the equilibrium function D = 
tanh(eVo/2fcBT), in agreement with the FDT. If the load 
conductance is large, Gi > Gm, the equilibrium noise 
of the load resistor may dominate in the circuit noise. 
In this case the function D may retain its equilibrium 
form, even if the mesoscopic conductor is biased. It is 
therefore interesting to consider the strong bias regime, 
eVia > kBT/R, so that the equilibrium noise contribu- 
tion can be neglected. Three important cases, which de- 
serve special consideration, are discussed below. 



o 
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I / / s 


\ rs 
1 r/ 




tanh{e/2Q) 




/ Z G M R=2 
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/ / 1 
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-**"""""""""""* s / / / 







e/Q 



FIG. 4: The symmetric detector function D(e, 0) is plotted 
versus normalized level spacing e/Q. for different values of the 
parameter TGmR- Typically, D is concave function of e, al- 
though for a super-Poissonian noise T > 1 it may become 
convex. Note that in the limit TGmR <S 1 the detector 
function has a power law behavior as compared to the expo- 
nential behavior of the equilibrium D = tanh(e/2f2) shown 
by the dashed line. 



A. Symmetric detector, A = 

This case is most relevant for the TP detector, which 
is symmetric detector. Using the zero-temperature limit 
F(uo) = u> in Eq. (29) and substituting the result to the 
Eq. (28), we obtain an important result 



for lei < eV M , 



TG M R{eV M -\e\) + \e\ 

(31) 

and D — ±1, otherwise. Here T = ^Z n T n {l — 
T n )/J2 n is the Fano factor of the system noise. Note 
that the slope of D at s = is equal to 1/ ' (bTGmRVm) = 
1/(217), where £1, we remind, is the circuit noise temper- 
ature. Interestingly, as we show below, this slope is uni- 
versal and same for an arbitrary mesoscopic conductor 
and arbitrary circuit. 

In the case of a very low load impedance, GmR "C 1, 
the result (31) can be written in the dimensionless form 



D(e,0) = 



e/n 



| e /n|' 



(32) 



It is plotted on Fig. 4 together with the equilibrium 
D(e, 0) = tanh(e/2f2). The non-equilibrium D has a 
power-low asymptotic at |e| — > oo, while the equilib- 
rium one shows an exponential behavior. Note also that 
for TGmR > 1 the detector function D(e, 0) is con- 
vex, which could be considered a signature of a super- 
Poissonian noise. 



8 



B. Asymmetric detector 

We consider a circuit far away from equilibrium, 
cVm > k B T/R. We expect that it might be difficult 
to adjust the DD detector precisely to cancel the asym- 
metry, therefore we first assume that the asymmetry is 
strong, |A| 3> cVm- Then, looking at the result (28), we 
expect that D = 1, i.e. the noise of the DD detector is 
Poissonian. In fact, more careful analysis shows that the 
strong asymmetry simply shifts the zero of the function 
D. This is because small or large value of the exponential 
e A/k B T ma y conl p ensa ted by the opposite effect in S 
due to the activation processes. 

Looking at the result (29) and (30) we find that 
for positive u) the dominant contribution is Sm(~^) = 
J r G M (^-eV M )e^ VM ~^ /kBT due to such activation pro- 
cesses. The contribution of the load resistor, Sl(—u) is 
small. On the other hand, both the load resistor and the 
mesoscopic system contribute to the term in Eq. (28): 
S(ui) = 2ui/R. Neglecting terms eVu and eVo compared 
to A, we finally obtain: 



D(e, A) = tanh 



eV D + sgn(A)^ M 



2k B T 



where the energy shift Em is given by 

E M = eV M + k B T\og{FG M R/2). 



(33) 



(34) 



Thus we arrive at the remarkable result that the only role 
of A is to fix the sign of the energy shift in (33). This 
fact is easily understood when we notice that the energy 
shift can be viewed as a drag or noise rectification effect, 
the direction of which depends on the sign of A. 

Interestingly, there is an additional contribution to Em 
in form of the logarithm, which contains the system Fano 
factor. It is exactly same parameter that also appears 
in the symmetric case (31). This additional shift may 
be interpreted as originating from high energy excita- 
tions that create the shot noise in mesoscopic system. Its 
explicit form depends on the assumption we made that 
the system is a coherent mesoscopic conductor. There- 
fore, it would be interesting to consider other examples 
of mesoscopic systems which may change the result (33) 
and (34). 



C. Classical noise regime 

So far we have discussed essentially quantum regime of 
the noise detection, where specific form of the function 
D depends on the choice of the system. In the rest of the 
paper we concentrate on the classical Markovian limit, 
which corresponds to a small detector bandwidth e, and 
demonstrate a number of universalities. 

We note that although the function D in Eqs. (31) 
and (32) behaves regularly at e = 0, it has been ob- 
tained in the limit e > IZfl using the result (16). If the 
detector bandwidth e tends to zero, P(e) given by Eq. 



(16) diverges and has to be replaced with the resummed 
version (21). Natural question arises is whether a consid- 
erable change in P(e), including appearance of the peak 
at e — 0, affects the symmetric detector function (31). 
The answer is no. We first check this for a fast and slow 
circuit limit, and prove the universality in the next sec- 
tion. 

Indeed, substituting either the function (21) or the 
function (23) into Eq. (24), for A = we obtain 



D(£,0) 



P(e) - P(s) 
P{e) + P(-e) 



e 

20' 



ife«0, (35) 



where, we remind, O is the circuit noise temperature (18). 
This result agrees with (31) as e — > 0. We stress however, 
that the result (35) is more general, since its derivation 
does not relay on the scattering theory 45 for a mesoscopic 
coherent conductor. 

We now in the position to investigate the effect of 
asymmetry. Restricting ourselves to the classical regime, 
e <C O, we write that generally P(e) = Po(e)(l + s/2Q,), 
where Po{e) — Po(—e) is the classical contribution. Sub- 
stituting this expression to the Eq. (24) , we arrive at 



D{e, A) = tanh(<5) 



20 cosh 2 (5)' 



2k B T 



• (36) 



Note that this result does not contradict the strongly 
asymmetric case (33), because here we assume £<SJ. It 
implies that at small asymmetry the zero of the detector 
function is shifted: 



D(e,A) = 



Eo 



20 



E 2 = (n/k B T)A. 



(37) 



Again, this shift is solely due to the second cumulant of 
current noise, and can be interpreted as an asymmetry 
induced noise rectification effect. This fact is important 
for the discussion in Sec. VII. 



VI. UNIVERSALITY OF CLASSICAL NOISE 
DETECTION REGIME 

In order to arrive at the result (35), we used Eqs. (21) 
and (23) for a fast and slow circuit, respectively. In gen- 
eral, the shape of the peak in P(E) depends on the circuit 
response time tq and circuit details via the impedance 
function Z(u>). In this section we show that, surprisingly, 
in the classical limit £«!! the detector function retains 
its form (35) and (36). 

We return now to the equations (12) and (14), and 
assume that the circuit impedance Z(u>) is arbitrary. The 
only requirement that we impose is that the interaction 
is weak: 1Z = GqZ{Q) <C 1. We focus on the long time 
limit of J2(t), so that the integral in (14) comes from 
small frequencies ui ~ TZfl, where the noise power can 
be approximately expanded as S(ui) = S(0) + d u S(0)u>. 
Consequently, J2(t) acquires two contributions that can 
be written as 

Mt) = G o S(0)H(t) + iG d u S(0)d t H{t), (38) 
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where 



H(t) 



duj 



\Z(u)\ 2 [l -cos(wi)]. 



(39) 



We are interested in time scales t ~ l/(TZCl), where 
the first term in (38) is of order 1, and the peak of the 
function P(E) is formed. Then the second term in (38) 
is of the order of 1Z, i.e. it is always small. Therefore, 
its contribution to the exponential in Eq. (12) should 
be expanded, giving the odd part of the P(E) function. 
Thus we obtain the following result: 

P(E) + P(-E) = it' 1 J dtexp[-G S(0)H(t)\ 

x cos(-Ei), (40a) 

P(E)-P(-E) = n^J dtcxp[-G S(p)H(t)] 

x G d u S(0)dtH(t)san(Et). (40b) 



It is easy to see that by the integration by parts the 
Eq. (40b) can be presented in the same form as Eq. (40a) . 
Thereby, independently of the exact function Z(cj), we 
arrive at the most general result for the classical noise 
regime, £<Cfi: 



P(E) - P(-E) _ d^SjO) 
P(E) + P(-E) ~ 5(0) 



E. 



(41) 



Using again the result (20), we arrive at the Eq. (35), 
which therefore holds for an arbitrary circuit. 



VII. THIRD CUMULANT CONTRIBUTION 

We have shown in the Sec. II that in the long time limit 
the quantum noise contribution to the correlator J2 (t) is 
small. The same remains true for the third cumulant, 
J 3 (t). Since the third cumulant contribution is small by 
the parameter 7t, we right from the beginning focus on 
its classical part, and rewrite Eq. (10) as follows: 



J 3 W = (1/6)([0(U-0(O)] 3 >. 



(42) 



Thus we see that J 3 (—t) = — J 3 (t). This breaks the sym- 
metry between right and left lead, Plr 7^ Prl, and the 
third cumulant adds to the potential difference across 
the tunnel junction. In the classical limit E <C O, where 
P(E) has a peak, this additional potential simply shifts 
the energy by a small amount E 3 that depends on the 
third cumulant of current: Plr(E) = P(E — E 3 ) and 
Plr{E) = P(E + E 3 ). Therefore, the function of the 
symmetric detector (35) has to be replaced with 



D(e,0) 



Plr{s) - P RL {-£) = e-E 3 
Plr{s) + Prl(s) 2fl 



(43) 



The same shift obviously takes place in the asymmet- 
ric case. However, there the shift E 3 adds to the shift 



E2 due to the noise rectification effect (sec the discussion 
in the end of Sec. VC). Fortunately, in contrast to the 
rectification shift, the energy E 3 depends on the direc- 
tion of current in a mesoscopic system. Therefore, ex- 
perimentally the third cumulant contribution can be ex- 
tracted by changing the direction of the current through 
the mesoscopic system. This experimental technique has 
been recently used to measure the third cumulant with 
the help of Josephson junction threshold detectors. 33 ' 34 
In addition, an important role in this context plays the 
universality of the Gaussian noise effect on the detector 
function D(e,0) proven in Sec. VI. In what follows we 
evaluate the shift E 3 for the cases of fast and slow circuit, 
depending on the circuit response time tc- 



A. Fast circuit 

In the case of fast circuit, TZflrc <C 1, the potential 
fluctuations are Markovian, so that Eq. (42) gives 



J 3 (t) = (e 3 /6)«V 3 K 



(44) 



where ((U 3 )) is the Markovian cumulant of the poten- 
tial. According to Refs. [10,14,15] it is given by ((U 3 )) = 
R 3 ((I 3 }}, where the total third cumulant of the current is 
equal to 

((/ 3 )) = ((i 3 M )) + endvSio) + \2{n/Rfd v R. (45) 

Here ((iff)) is the intrinsic third cumulant of the system 
current, and the second and third terms are the "environ- 
mental" and nonlinear cascade corrections, respectively. 
They originate from the circuit backaction. 

It is useful to write J 3 in the form that explicitly shows 
the coupling constant: J 3 (t) = (l/6)(2nK/e) 3 {(I 3 )) t. We 
see that indeed such a contribution to the correlator sim- 
ply shifts the energy in the Fourier transform (11) for the 
probability distribution functions by the amount 



E 3 



(l/6)(27rft/e) 3 ((/ 3 )), ftfiT C «l. 



(46) 



In order to estimate the relative effect of the third cu- 
mulant we note that the width of the peak in P(E), 
where the detector signal is maximum, is of the order 
of 1Z£1, so that D <~ 1Z. The energy shift can be esti- 
mated as (TZ/e) 3 ((I 3 )) - TZ 2 G M R^- Therefore, the rela- 
tive contribution of the third cumulant is of the order of 
KG m /{Gm + G l ) « 1. 



B. Slow circuit 

In the case of slow circuit, TZflrc ^> 1, the detector 
"feels" slow fluctuations of the potential. Therefore, one 
can approximate <p(t) — </>(0) = edVt, according to exact 
calculations in the Sec. IIIB. Then the equation (42) 
gives: 



J 3 (t) = (e 3 /6)((SVf)t 3 



(47) 
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where ((SV) 3 } is the third cumulant of equal time fluc- 
tuations of the potential V. In Refs. [12,13] it has been 
shown that ((5V) 3 ) = (R/3C 2 )((I 3 )}, where the total cur- 
rent cumulant in this case is given by 

((i 3 )) = + md v s(o) + 3(n/R) 2 d v R. (48) 

It contains the intrinsic cumulant of the system current, 
((/£f)), and the cascade corrections. Note that in this 
case, IZQtc 3> 1, the cascade corrections are smaller 
compared the those for a fast circuit, see Eq. (45). This 
fact has been recently experimentally verified in Ref. [12]. 

We now substitute the small term (47) to the definition 
(11). This gives P LR = [1 + (e 3 /6)((5V) 3 )d 3 E ]P(E) and 
Prl = [1 - {e 3 /6)({SV) 3 )d 3 E }P(E). Using Eqs. (23) and 
(43) we obtain 

^=!S5« /3 »< ™ tc>>1 - (49) 

Note that this energy shift is smaller than the one for 
the case of fast circuit, Eq. (46), by the parameter 
1/(R£It c ) < 1. Since the width of the distribution P{E) 
is of the order of ^TZCI/tc, the relative contribution of 
the third cumulant to the function D can be estimated 
as KG M /[(G M + G L )(KnT C ) 1/2 } « 1. 

VIII. OUTLOOK 

We have presented the theory of quantum noise detec- 
tors based on the resonant tunneling phenomenon. It is 
summarized in the introduction, which can also be used 
as a guide to most important results. Here we briefly 
discuss related problems which yet to be solved. First 
of all, it would be interesting to relax the condition of a 
weak coupling. In the case 1Z <~ 1 fluctuation dissipation 
relations may contain an information about the full dis- 
tribution of the fluctuating potential. Interestingly, it has 
been shown in Ref. [46] that the double dot system in the 
adiabatic regime, tc(Id) 1, may serve as a nonlinear 
element which generates an instability in the mesoscopic 
circuit. It then may be used as on-chip threshold de- 
tector of rare event in transport. The difficulty of this 
problem is that the dynamical Coulomb blockade effect 
in this case is not generally negligible. 

In the case of a TP detector the excitation of an 
electron-hole pair in the QPC may cause a transition in 



the DD system. This competing quasiparticle process re- 
duces the precision of the detection of collective charge 
excitations in the measurement circuit. Intuitively, one 
should keep the current through the QPC on a very low 
level. However, this will reduce the rate of the measure- 
ment. Moreover, the quasiparticle process is interesting 
in itself and should be investigated theoretically. 

We think that the physics of double-dot systems de- 
scribed here is rather universal and should be same in 
various two-level systems of a different nature. Never- 
theless, it is important to consider other systems too. 
Moreover, it would be interesting to generalize present 
results to the case of a quantum detector with many lev- 
els with the energies e n , n — 1,2,.... There is a hope 
that such system will be able to detect high order corre- 
lators of current at finite frequencies equal to the energies 

Concerning specific results presented in the paper, two 
problems remain to be solved. First, we have shown that 
in the classical noise detection regime the effect of the 
second cumulant of the system current is universal, i.e. it 
does not depend on the circuit details. On the contrary, 
the third cumulant contribution depends on the circuit 
response time and has been found here in the limit of 
fast and slow circuit. Interesting problem, which may 
also be experimentally very relevant, is to find the third 
cumulant contribution including cascade correction for 
arbitrary circuit. 

Second, the most dramatic effect of a non-equilibrium 
system noise on the detector function is that the expo- 
nential behavior (1) is replaced with the power law func- 
tions (31) and (32). Thus the power law behavior is a sig- 
nature of non-equilibrium processes. However, this our 
result has been obtained by considering a coherent non- 
interacting mesoscopic conductor as an example of the 
system. Therefore, it would be interesting to consider 
other systems in order to check the generality of our con- 
clusion. 
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